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Abstract 

The complex nature of inertial confinement fusion (ICF) experiments results in a very large number of experimental 
parameters that are only known with limited reliability. These parameters, combined with the myriad physical models 
that govern target evolution, make the reliable extraction of physics from experimental campaigns very difficult. We 
develop an inference method that allows all important experimental parameters, and previous knowledge, to be taken 
into account when investigating underlying microphysics models. The result is framed as a modified analysis which 
is easy to implement in existing analyses, and quite portable. We present a first application to a recent convergent 
ablator experiment performed at the NIF, and investigate the effect of variations in all physical dimensions of the 
target (very difficult to do using other methods) . We show that for well characterised targets in which dimensions vary 
at the 0.5% level there is little effect, but 3% variations change the results of inferences dramatically. Our Bayesian 
method allows particular inference results to be associated with prior errors in microphysics models; in our example, 
tuning the carbon opacity to match experimental data (i.e., ignoring prior knowledge) is equivalent to an assumed prior 
error of 400% in the tabop opacity tables. This large error is unreasonable, underlining the importance of including 
prior knowledge in the analysis of these experiments. 

Keywords: inertial confinement fusion, radiation hydrodynamic simulation, Bayesian inference, plasma opacity, 
uncertainty analysis, convergent ablator, national ignition facility 



1. Introduction 

The design of experimental schemes to reach thermonu- 
clear ignition and burn in laser driven targets involves 
complex models that incorporate many physical effects 
[1] . The radiation- hydrodynamic simulations used to pre- 
dict the evolution of fusion capsules [2] therefore contain 
a huge number of physical parameters which play an im- 
portant role. The resulting laser targets are correspond- 
ingly complex, with a large number of design parameters. 
In a typical inertial confinement fusion (ICF) experiment 
performed at large laser facilities such as the national ig- 
nition facility (NIF) [3], there are many tens of variables 
that play an important role in determining target evolu- 
tion [3]. This poses a difficult problem for data analysis 
since these parameters should not be neglected but are too 
numerous to treat directly using the standard methods of, 
for example, particle physics where Monte Carlo sampling 
of noise sources is often used [S] . In this paper we develop 
a method that allows all important variables to be in- 
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eluded, along with prior work on microphysics models, in 
a consistent and efficient analysis. The approach has been 
designed to couple with existing radiation-hydrodynamics 
simulation codes without modification; in fact simulations 
are treated as a 'black box' making the method applicable 
to a large class of difficult data analysis problems. This 
approach also allows us to avoid the complex fitting func- 
tions used in other approaches [HI [7j , which are unlikely 
to capture the complex behavior of ICF experiments close 
to ignition (and are unsuitable for such large problems in 
any case). 

The data analysis approach that we describe is particu- 
larly important when considering the results of the recent 
national ignition campaign (NIC). Throughout the cam- 
paign, post-shot simulations failed to match the observed 
data; the implication is that simulations, or their under- 
lying microphysics models, are inaccurate. Determining 
which of the models should be investigated, and produc- 

38 ing a consistent picture of the implied error, is a difficult 

39 task and forms a major motivation for this work. 

40 In fact, modifications to various physical parameters, 

41 even unrealistically large modifications, often cannot pro- 

42 duce a match to NIC data. In this situation the neglect 

43 of important variables and prior knowledge has a dra- 

44 matic effect on the results of inference (even if they are 
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45 purely sources of noise). Our method allows these effects 

46 to be included with almost no computational overhead. 

47 We demonstrate this by presenting an analysis of a single 

48 NIC convergent ablator (conA) shot [HE], N110625. We 

49 find that variations in the dimensions of the target can 

50 have a dramatic effect on the inferred drive and carbon 

51 opacity, although this is mitigated by thorough metrology 

52 of the target. Our method also allows prior knowledge to 

53 be included and in the case of ICF we find that this is an 

54 extremely important factor. We use it to investigate the 

55 implied error in microphysics models associated with ne- 

56 glecting this prior work. The inclusion of this prior knowl- 

57 edge is an important strength of the Bayesian method, as 

58 it provides context for observed data and therefore allows 

59 meaningful information to be inferred even from a sin- 

60 gle experimental result. The total information from a set 

61 of experiments can be viewed as a series of such single- 

62 shot inferences, allowing the analysis performed here to 

63 be generalised to full experimental campaigns very easily. 

64 The approach is to treat the output of the simula- 

65 tion code as probabilistic, and to apply standard meth- 

66 ods of Bayesian analysis |10| . The probabilistic nature 

67 of simulations is due to variations in the myriad inipor- 

68 tant variables (or 'nuisance parameters'). We derive a 

69 semi-analytic expression in which the dependence on in- 

70 teresting physics is retained but all other variables are 

71 represented by an analytic information loss. The result 

72 is framed as a modified analysis which is easy to im- 

73 plement, portable, and allows all available data to be in- 

74 eluded in a single analysis. 

75 We begin by elaborating on the challenges we have 

76 already introduced. We then develop our inference ap- 

77 proach in section [3j and discuss methods for its applica- 

78 tion in section [4] Finally, the importance of including 

79 all important variables and prior knowledge is demon- 

80 strated with an example application to a single NIC shot, 

81 N110625. 

82 2. Challenges for data analysis from ICF experi- 

83 ments 

84 In current analyses, particular data (chosen largely 

85 through experience) are preferentially matched by vary- 

86 ing inputs that are considered to be unreliable, such as 

87 X ray drive [TTJ [T^] . This approach has been very useful 

88 in testing the consistency between simulations and ex- 

89 periment, however it is potentially sensitive to noise and 

90 gives little information about the physical origin of incon- 

91 sistencies. Increasing the number of inferred parameters 

92 is essential to gaining more information about underlying 

93 physics models. 

94 Radiation-hydrodynamic simulations represent a non- 
95 linear map from the space of physical models that we 

96 wish to investigate to the data that are collected in an 

97 experiment. The nature of the simulations often means 

98 that they are not amenable to adjoint differentiation [13,, 



99 are discontinuous, and may be noisy; these complex fea- 

100 tures can make standard methods of searching the space 

101 of physical parameters quite unreliable. This limits the 

102 number of parameters that can be reliably inferred. The 

103 nuisance parameters included by our method result in a 

104 smoothing of the code output, allowing the use of ad- 

105 vanced methods and an increase in the number of physical 

106 parameters that can be investigated. 

107 We have already described the difficulties associated 

108 with the large numbers of target parameters involved 

109 in ICF experiments. Although many of these are con- 
no strained by manufacturing precision and target metrol- 

111 ogy, it has already been seen that their large number can 

112 have an important impact on the output of simulation 

113 codes [4j . There will be a corresponding effect on infer- 

114 ence results, and we aim to investigate this. 

115 The physical parameters we aim to infer often refer to 

116 microscopic physics (for example opacities or equations 

117 of state) that are understood using other, separate, com- 

118 puter simulations. These simulations are highly complex 

119 and have been investigated both theoretically and experi- 

120 mentally for many decades; the expected systematic error 

121 bars on their outputs are therefore quite small. This error 

122 bar plays an important role in data analysis by ensuring 

123 that the results are physically reasonable, and this moti- 

124 vates our Bayesian approach. 

125 3. Probabilistic output from a deterministic sim- 

126 ulation code the importance of nuisance pa- 

127 rameters 

128 The fundamental problem is to develop a method of 

129 exploring the huge space of parameters that can affect 

130 the outcome of a simulation. As discussed, in the case 

131 of ICF data there is no point in this space for which all 

132 data are correctly simulated. In general there may be a 

133 set of points that give comparable agreement. The best 

134 fit is found by defining a figure of merit that takes into 

135 account the difference between observed and simulated 

136 values of all data points, as well as the difference between 

137 simulation parameters and the expected physical reality. 

138 In this section we outline a figure of merit that is based 

139 on the Bayesian posterior probability of a point in phase 

140 space (the maximum a posteriori^ or MAP, solution |10|). 

141 and use an analytic prior-predictive approach to reduce 

142 the phase space to manageable size. 

143 We begin by splitting the set of all parameters into two; 

144 • 'Interesting Parameters' 9 - Physically significant pa- 

145 rameters that we aim to infer from experiment data. 

146 For example material equation of state, opacities, 

147 conductivities, 

148 • 'Nuisance Parameters' rj - Other parameters that 

149 have an effect on simulations but are not of direct 

150 physical significance. These are usually known with 
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good precision, for example target dimensions, laser 173 of nuisance parameters results in a ~ 20 dimensional inte- 
powers, 174 gral, this is prohibitively expensive. Even if the radiation- 

175 hydrodynamics can be modelled by some fast surrogate 
In our model, inference is performed on the interesting ^^^^j ^ Gaussian Process or through other techniques 
parameters only Bayes' theorem allows us to write down [JT] ) , which itself is very difficult given the size of 

the probability distribution of the interesting parameters ^^e space we must consider, the integral is still too ex- 
once the experiment has been performed (the posterior), pensive. We instead evaluate the integral by assuming a 
in terms of the probability distribution before the experi- ^j^^^^^. response to nuisance parameters, 
ment (the prior) and the probability of the experimental 

data (the likelihood). Bayes' theorem is dm{0,r]) — djn{9,ri = r/o) + A{r]Q — 77) . (2) 

P{dexp\d)P{0) In the above 770 is the nominal value of the nuisance pa- 

P{S\dexp) — P[dex ) ^'^^ rameters, typically zero. The linear response matrix A 

ffd dd P(d d 9 ) las can be populated using a small number of simulations; 

= — — — — — — ^ , 184 once this has been done, the matrix A is entirely portable 

\"'<i^P' 185 and may be used in all subsequent analyses of this type 

where de.p is the vector of experimental data and we ^i*^"'^* ^^^^l^^^' calculation. 

have introduced the code output d„ and the nuisance ^he case of linear response, equation is very useful 
parameters as marginalised variables. This allows us to ^^l^^^ ^^^^1^*^^ treatment of the nuisance parame- 

introduce the known measurement error and prior distri- t^^^^- Assuming that the experimental measurement errors 
butions of the nuisance parameters later. Such an ap- "° nuisance parameter variations are described by un- 

1 • • 1 J. i • 4.U 4- • 4. 1 J J. 191 correlated normal distributions with correlation matrices 

proach is equivalent to assuming that experimental data 

are the simulation results plus a randomly distributed '^^p ^"^^ 

error, as is done in other approaches [Hill]. Writing r g-(de.p-d„(e,,)))^A-Jj,(de.p-d„(e,^)) 

P{dexp,dm,9,r]) = P{dexp\dm,d,r])P{dm\9,ri)P{9,r]) and P{dexp\9) ^ dr] I ia ^F 

introducing the deterministic nature of the simulation vi^"'/ "^''l-^expl 

code, p-iv-vo)"^ iv-va) ' 

P{d^\9,7^)=5{d^-d^{9,i^)) , > 



X 



V(2^)"''|A,| 



the integration over d,n can be performed trivially. The 
result is 



the result is 



_ P{^) i^-j ^--j^-j In the above, d,n{9) = d„i{9,r]f)) is the simulation result 

P{dexp) "^^^ ' for nominal nuisance parameters and the matrices a and 
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f3 satisfy the equations 

,T aT, 



The likelihood P(c?e2;pMm(^i implicitly contains the 

experimental error distribution and the code output as a a a — A A^^pA + A,^ 

function of all parameters dm{9, nf). The two components a = A 

of the prior distribution P{9, rj) = P{9)P{rj) describe the 



157 expected distributions of nuisance and interesting param- im These expressions are the multivariate generalisation of 

158 eters before the experiment has been performed; these 195 the usual quadrature error propagation formula; it should 

159 are determined by the experimental design, target man- i96 be noted that even if nuisance parameters and experimen- 

160 ufacturing tolerances, previous experimental results and 197 tal errors are independent to begin with (i.e., if Ae^p and 

161 expert opinion. 198 A^ are diagonal) , the response of the simulations means 

162 Equation ([!]) describes the relationship between the 199 that the likelihood can become strongly correlated. These 

163 probability distributions of the interesting parameters be- 200 potentially strong correlations arise due to the determin- 

164 fore and after an experiment. The details of the relation- 201 istic nature of the simulation code and play a very impor- 
ship are approximated by the simulation code, and con- 202 tant role in the inference procedure described in the next 
tained in the likelihood function. Data analysis, then, is 203 section, 
based on the evaluation of the integral in the definition 

of the likelihood. In its general form this involves the Inference of interesting parameters from exper- 

169 integration of simulation output over the entire nuisance imental data 

170 parameter space; it is common to evaluate this using a 

171 Monte-Carlo sampling of the integrand (see, for example, 206 The results of the previous section allow the efficient 

172 [161 ). For our application, where even a conservative set 207 calculation of the likelihood as a function of interesting 
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208 parameters, without neglecting other important variables 256 In the following section we present an example applica- 

209 or prior knowledge. As discussed in section [2j this can be 257 tion of the method. For simplicity the parameter space is 

210 expected to give a significant improvement in data anal- 25s small allowing the likelihood and posterior to be explored 

211 ysis results. The marginalisation of nuisance parameters 259 directly. The advanced techniques discussed above are 

212 represents an averaging that smooths the response of sim- 260 therefore not needed. In a forthcoming paper we con- 

213 ulations, making them more well-behaved. This allows us 261 sider more complex cases for which we develop a genetic 

214 to use standard numerical techniques. 262 algorithm that is optimised for the sparse datasets en- 

215 The best fit to data, taking into account nuisance pa- 263 countered in ICF research. 

216 rameters and prior knowledge, is given by the parameters 

21. that maximise the posterior probabihty P{e\d,,p) (see ^ Application to NIF convergent ablator experi- 

218 equation It is convenient to minimise the informa- mental data 

219 tion, I{9\dexp) = —LogP(9\dexp), which using equations 

220 and (|3| is 266 In order to demonstrate the application to an actual 

^ 2 267 inference problem, we now consider experimental data 

l{0\dexp) — ^L_5^PiL_ — 268 taken on a NIC convergent ablator (conA) experiment 
i ^exp,i 269 [5J E . In this design, an ICF capsule is imploded and 

- (4xp - d„,{9)f Pidexp - d„,i9)) =™ backht by emission from a nearby high Z plasma. This al- 
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271 lows time- and space- resolved measurement of the plasma 
, ^-In (|A^||a^Q!|) — lnP(6') . (4) 272 density during the implosion. This is analysed to give 

273 time resolved measurement of fuel shell position, velocity, 

221 The above equation has the form of a modified func- 274 and line density. Simple models show that these quanti- 

222 tion, and is derived by assuming that A^xp is diagonal. 275 ties are sensitive to the details of the X ray drive from 

223 Note that the dependence of the first term on 9 through 276 the hohlraum, and to radiation transport in the capsule 

224 the simulation dm{9) means that even in the absence of 277 ablator ^j. 

225 nuisance parameters and prior knowledge the likelihood is 278 In these experiments the measured implosion velocities 

226 non-normal. Equation Q can be interpreted as an infor- 279 are consistently lower than simulated predictions, possi- 

227 mation processing rule jj] ; the first 3 terms on the right 280 bly suggesting a reduced X ray drive. Absorption of the 

228 hand side are the information gained from the experiment, 281 drive X rays by carbon in the ablator plastic also plays an 

229 and the final term is the information about the interesting 282 important role, and simulations show that an increased 

230 parameters before the experiment was performed. In that 283 carbon opacity can improve agreement [Mj. These pa- 

231 sense it is clear that the positive definite matrix rep- 284 rameters can be used to tune simulations to agree with 

232 resents a loss of information due to nuisance parameters. 285 experiment, however such an approach runs the risk of de- 

233 As mentioned, once /3^/3 has been computed the evalua- 286 stroying the predictive capabilities of codes when run far 

234 tion of the modified oi^ly requires a single simulation. 287 from existing experimental data (a common occurrence in 

235 In an actual inference problem we are interested in the 288 all areas of HEDP, and a possible problem when design- 

236 values of 9 that give the best fit to the experimental data. 289 ing improved ICF targets). We aim to analyse the signif- 

237 This requires the numerical minimisation of equation Q . 290 icance of these experimental results with respect to the 

238 The well behaved nature of the marginalised likelihood al- 291 drive and carbon opacity by inferring the values of modi- 

239 lows us to use standard methods; two common approaches 292 fiers to those quantities in the presence of many nuisance 

240 are 293 parameters and of prior knowledge. The prior distribu- 

241 . Markov Chain Monte Carlo (MCMC) - This ap- *^°^f P^^^^ ^^^'^l multipliers are interpreted 

1 . • i- i ii i- , 295 as the uncertainties in the on-line calculations of opacity 

242 proach gives an approximation to the entire poste- d d ' 

243 rior information [19]. This is extremely useful to mi • • 1 i 1 tt^tt^t^ « t • 

,1 1 .• c • -c J , 297 ihe inference is based on the HYDRA radiation- 

244 the evaluation of error bars on interred parameters. . __ 

m, , J rr ■ ii i ii ii J • i ■ 298 hydrodynamics code 20 . The parameters of interest are 

245 1 he trade-ort IS that these methods require extensive / ' — , . ,. 

,1 . , . J J j-o: li i • 11 1 299 the values of two dimensionless multipliers; one is applied 

246 burn m periods and are ditncult to run m parallel. / iT • 1 1 

T 1- 1 ■ 1 • 1 i- 300 to the X ray drive spectrum (at all times and photon 

247 in our applications where a single simulation repre- • \ 1 i 1 . V 1 1 

.J, , ... , , 1 ii • • 301 energies) and the other is applied to the carbon opac- 

248 sents a signihcant computational overhead, this is a . ' i • ■ 1 i . n 

J. J , 302 ity (all temperatures, densities and photon energies). We 

249 major disadvantage; , . . n . n 

303 take into account 29 nuisance parameters, allowing all 

250 • Genetic Algorithm (GA) - This method uses ideas 304 capsule dimensions, material densities and material com- 

251 taken from genetics to efficiently find the minimum 305 positions 4 to vary. These parameters are allowed to 

252 of a function. It is very easy to run in parallel and 306 vary according to two distributions with standard devia- 

253 so is well suited to our application. The final result 307 tions of 0.5% and 3% respectively. These represent well 

254 is the position of the minimum only and so some 308 constrained target parameters (many NIC capsule dimen- 

255 approximation is required to calculate error bars [10 , 309 sions are known to better than 0.5%), and ones with 
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310 more uncertainty (demonstrating the potential effect of 

311 nuisance parameters on inferred physics). For this rel- 

312 atively small inference it is feasible to generate a set of 

313 simulations that span the 2D parameter space. For each 

314 point, defined by the multipliers {Aj_rive, '-^Copacjj we run 

315 HYDRA and extract the implosion velocity, ablator mass 

316 fraction, and time at which the implosion reaches a radius 

317 of 310/im. These quantities are compared to experimental 

318 values taken from radiography |21| . 

319 In figure [l] we plot the information in the likelihood as 

320 a function of ^drive and Acopac, calculated using equa- 

321 tion ([3]) with different values of the modification matrix 

322 /3"^/3. These plots represent our modified when the 

323 prior distribution P{0) is neglected. In (a) nuisance pa- 

324 rameters are neglected (/3-^/3 = 0), and in (b) the modi- 

325 fication is calculated as described for all 29 nuisance pa- 

326 rameters varying at the 3% level. In the case with 0.5% 

327 variations, nuisance parameters have a small effect and 

328 the likelihood is very similar to figure [ij a). The positions 

329 of the minima are marked with red points, making the ef- 

330 feet of nuisance parameters clear. This shift in minimum 

331 is very important in the subsequent analysis. 

332 To further quantify the differences we perform a set 

333 of inferences based on the calculated likelihoods. The 

334 specific choice of prior distribution is often a difficult is- 

335 sue since it can be a subjective choice that has a direct 

336 influence on inference results. For this reason we per- 

337 form a range of inferences with prior distributions for 

338 {Adrive, Acopac) of Varying width. This allows us to take 

339 into account the dependence of the inference results on 

340 the prior, and place limits on the actual prior for various 

341 results. 

342 We begin with a reasonable estimate of the uncertain- 

343 ties in opacity and drive models of 10% and 20% respec- 

344 tively. This defines our nominal prior as a normal distri- 

345 bution centered on (1, 1), with covariance matrix 



0.1^ 




0.22 



(5) 



346 A set of inference results are found by scaling this covari- 

347 ance, thereby changing the assumed prior error in micro- 

348 physics models (and the relative importance of prior and 

349 experimental information). For a very large scaling of ([5|, 

350 the prior is flat and our analysis reproduces the maximum 

351 likelihood (ML) result; for a small scaling factor the prior 

352 tends to a (5-function and the minimum of equation Q is 

353 at {Adrtve, Acopac) = (1, 1) (their prior values). 

354 In figure [2] we plot the trajectories of the best fit as 

355 the prior covariance is scaled from small to large. The 

356 trajectory for calculations that neglect nuisance parame- 

357 ters, and that include them at the 0.5% level, overlay each 

358 other and are shown in purple; note the slight shift in the 

359 ML result at the right hand end. The 3% case is plotted 

360 in green. The shapes of the trajectories are determined 

361 by all the factors we have discussed so far, not least the 
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(a) No nuisance parameters 
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(b) Linear model for nuisance parameters with 3% variations 

Figure 1: Information in the likelihood for multipliers 
placed on the carbon opacity and X ray drive for a NIC 
conA experiment. Panel (a) shows the result when no nui- 
sance parameters are included, and (b) shows the effect 
of including target metrology as nuisance parameters. 
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Figure 2: Trajectories of the best fit to experimental data 
from a NIC conA experiment, as the prior width is varied 
(see equation ([s])). The blue line shows the result when 
nuisance parameters are ignored, or included at the 0.5% 
level. The two red points at the right hand end represent 
the maxima of the likelihood for these two cases. The 
green line shows the case when nuisance parameters are 
included at the 3% level. As the prior is scaled from 
a (5-function, through our best estimated defined by ([5|, 
to flat, the inferred results tracks from the prior results 
(1, 1) to the minimum of the likelihood functions plotted 
in figure [T] The figure also shows contours that define a 
change in multiplier of 5% from each end point. 



362 shape of the likelihood (i.e., the effect of nuisance param- 

363 eters). The left hand end of the trajectories corresponds 

364 to small prior error and reproduces the prior result. The 

365 right hand end of each line is the flat prior result; as we 

366 have already seen in figure [l] the inclusion of nuisance pa- 

367 rameters at the 3% level has a very significant effect on 

368 the inferred values of our interesting parameters. 

369 The wide difference between the start and end points of 

370 all trajectories in figure [2] clearly shows that the prior dis- 

371 tribution has an extremely important role to play in our 

372 analysis. For our nominal prior, defined by the covariance 

373 ([5]), we find that the prior is in fact more important than 

374 the details of the nuisance parameters regardless of their 

375 distribution widths, giving inference results that are al- 

376 most the same; {Adrive, ^Copac) = (1.03,0.94). The ML 

377 analysis, that neglects the prior, will then result in a sig- 

378 nificantly different result. This is true even for extremely 

379 broad priors; for our MAP analysis (which includes both 

380 prior and nuisance parameters) to reproduce the 0.5% nui- 

381 sance parameter ML result to within 5% (shown by the 

382 dashed contours in figure [2| , the prior covariance must 

383 be scaled so that the prior errors in opacity and drive are 

384 more than 400% and 800% respectively. The simulations 

385 on which the opacity and drive are based can be expected 

386 to be much more accurate that this, giving further sup- 

387 port to the importance of the prior. 



388 6. Discussion and Conclusions 

389 We have developed a Bayesian model for investigation 

390 of underlying physics using complex HED experiments. 

391 The model allows for the inclusion of complications aris- 

392 ing in experiments by using an approximate description 

393 of so-called nuisance parameters, and of previous investi- 

394 gations through a Bayesian prior. The result is a modi- 

395 fled function that can be easily incorporated into any 

396 analysis using standard methods. This approach allows 

397 complex simulations to be treated as black box transfor- 

398 mations from physical models to experimental data and 

399 so is suitable for application in a wide range of physi- 

400 cal applications. The linear response model described is 

401 the basis of the usual 'Normal Linear' model [10] • How- 

402 ever, unlike that model, the use of complex simulations 

403 to describe interesting parameters and the resultant cor- 

404 relations between nuisance parameters results in a non- 
405 normal posterior. 

406 In the case of ICF experiments, the linear response 

407 approximation may not be sufficient. The difficult task 

408 of achieving thermonuclear ignition requires that target 

409 designs are highly optimised; a change in nuisance pa- 

410 rameters in either direction is likely to produce a reduc- 

411 tion in target performance. Such nonlinear behavior can 

412 be important, and is not captured by the current ap- 

413 proach. Test calculations for a reduced problem, includ- 

414 ing quadratic response to nuisance parameters, suggests 

415 that these effects are significant in the analysis of ICF 

416 data. A major piece of further work is to develop an effi- 

417 cient way of including nonlinearity. 

418 In the final sections of this paper we have applied our 

419 analysis to a single NIC experiment. We attempt to de- 

420 scribe deficiencies in radiation transport physics through 

421 multipliers on two physical quantities, and infer the poste- 

422 rior values of these multipliers. This process is a common 

423 one in the analysis of NIF data, and is usually viewed as 

424 the tuning of simulations to allow more reliable target de- 

425 sign. In this work we interpret the results of this process 

426 as a measure of the uncertainty in the underlying physi- 

427 cal models, which are often applied in regimes where they 

428 are untested. Only by improvement of these models, mo- 

429 tivated by the kind of data analysis described here, can a 

430 truly predictive simulation be developed. 

431 The particular example given here is sufficient to 

432 demonstrate the importance of an integrated approach 

433 to data analysis, and provides compelling evidence that 

434 a straightforward fit to experimental data, ignoring prior 

435 knowledge, can give misleading results. For the very well 

436 characterised targets used at the NIF, certain dimensions 

437 are known to better than the 0.5% accuracy we allow in 

438 this work, however other nuisance parameters (for exam- 

439 pie material densities) could vary over a larger range. We 

440 have demonstrated that these nuisance parameters may 

441 have an important effect; our method allows a complete 

442 description of the problem. Alongside the nuisance pa- 
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rameters that we have included in this demonstration, 
there are also many other simulation inputs which can be 
treated as nuisance parameters in the same way. 

We demonstrated a novel method of analysing the im- 
portance of prior knowledge by referencing the possi- 
ble conclusions from data to limits on prior distribu- 
tion widths. The multipliers used here do not, however, 
provide an insight into specific problems in underlying 
physics; it is also true that these multipliers only describe 
the average modification to theory that is required. Any 
inferred physical modifier will lose its meaning when the 
simulations used in the inference have other unknown in- 
accuracies, and this is certainly the case in our first ap- 
plication. We begin addressing these problems in a forth- 
coming paper. 

The work presented here represents the first steps to 
providing a clearer view of problems with physics mod- 
els from experimental data, in cases where the experi- 
ments are very complex. Although we concentrate on 
ICF experiments here, nuisance parameters can be ex- 
pected to be important in all HED experiments, in partic- 
ular those where target plasmas are less well constrained. 
The portability of our method makes its application to 
other experiments very easy. The computational frame- 
work described also provides the opportunity for Bayesian 
experimental design [22\, allowing future experiments to 
provide a significant measurement of difficult aspects of 
underlying physics [23] ■ The integrated approach that 
we propose may also facilitate discovery of new rules and 
phenomenology that govern the evolution of these com- 
plex systems. 
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